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ABSTRACT 

At the conclusion of cryogenic vacuum testing of the James Webb Space Telescope Optical Telescope Element 
Integrated Science Instrument Module (JWST-OTIS) in NASA Johnson Space Center’s (JSCs) thermal vacuum (TV) 
Chamber A, contamination control (CC) engineers are postulating that chamber particulate material stirred up by the 
repressurization process may be kept from falling into the Integrated Science Instmment Module (ISIM) interior to some 
degree by activating instrument purge flows over some initial period before opening the chamber valves. This 
manuscript describes development of a series of models designed to describe this process. The models are strung 
together in tandem with a fictitious set of conditions to estimate overpressure evolution from which net outflow velocity 
behavior may be obtained. Creeping flow assumptions are then used to determine the maximum particle size that may 
be kept suspended above the ISIM aperture, keeping smaller particles from settling within the instrument module. 

Keywords: contamination control, purging, backfill, repressurization, venting 


1. INTRODUCTION 



Figure 1. Illustration of the James Webb Space Telescope. 1 

At the end of James Webb Space Telescope (JWST) OTIS (Optical Telescope Element -OTE-Integrated Science 
Instmment Module-ISIM) Cryovac testing in NASA Johnson Space Center’s (JSC’s) thermal vacuum (TV) Chamber A, 
contamination control (CC) engineers are mooting the idea that particulates stirred up by the repressurization process 
may be kept from reaching the ISIM interior to some degree by first activating instrument purge flows for a set period 
before opening the chamber pressurizing valves. The ISIM aperture faces upward along the TV chamber’s vertical axis. 
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This paper describes development of a series of models designed to describe the overall process as the gas transits 
various domains of rarefaction and other phases. These models are strung together in tandem to estimate overpressure 
evolution from which net outflow velocity behavior may be obtained. Creeping flow assumptions due to Stokes 2 are 
then used to determine the maximum particle size that may be lofted, keeping smaller particles from settling within the 
ISIM compartment. 


2. MASS CONSERVATION STATEMENT 



Figure 2. Schematic representation of a simple rigid vessel venting system. 

For this application, the rate of mass accumulation for a gaseous species within rigid ISIM volume V is made up of the 
rate of mass generation iii ln due to purging minus the net rate at which mass flux exits across that volume’s bounding 
surface S. The system is visualized in Fig. 2. Time is represented by t, mass density by p, and local flow velocity by U. 

~r I P dV = m m — § pUJ ■ dS . (1) 

dt v s 

Assuming thermodynamic conditions within this unit may be lumped, then the volume integral may be replaced by V. 
For an ideal gas under isothermal conditions, Eq. (1) may be simply recast in terms of pressure p, 

V — = /h jn RT-§pUdS. (2) 

dt 5 

In Eq. (2), species gas constant R for average molar mass M is related to the universal gas constant R u by R = R U /M. 
Considering molecular nitrogen as the purge gas, M = 28 g/mole. The equation now has units of gas load Q, often 
presented in Torr-L/s or W . 3 Since in = pG where G represents a local volumetric flow rate. 


Q = riiRT = pGRT = pG . 


( 3 ) 


2.1 Venting 


Where venting occurs across k orifices in parallel with areas A k against a background pressure p b , the surface integral 
may be replaced by a sum of K conductive terms. 3 


K 


V^- = m m RT-±F k {p-p b ). 
dr k 


(4) 


Conductance F has units identical to volumetric flow rate, but is defined using a pressure difference across the vessel 
boundary. For a single vent out the ISIM aperture, one may dispense with the summation notation. 

For states 1 & 2 describing conditions within V and the background outside, 3 

Fl-2 (5) 

Pl~P2 P\~P2 

Equation 5 will be used to develop relevant venting terms as the repressurization process transits different flow regimes. 


2.2 Venting Equations 

Analysis of this arrangement actually requires describing a number of stages tied together over time in series by various 
physical constraints. By purging through the instruments, this feature is meant to ensure the pressure will always be 
higher within the ISIM volume than outside of it even once the vacuum chamber valve is opened (an assumption 
reviewed in Section 3.4). A positive overpressure will produce a net outward flow across the upward-facing instrument 
module aperture, preventing particles below a certain size from dropping into ISIM. The average velocity associated 
with this overpressure may be determined by application of an energy conservation statement. If the net flow velocity is 
exceedingly low however, one may determine that the maximum particle size lofted against gravity will be so small that 
the purge application would be considered ineffective. 

In this application it is assumed that the instrument purge is initiated one hour before opening TV chamber 
repressurization valves. The aggregate purge flow rate G/ through the instruments is fixed and the ISIM volume is 

denoted by Vj . The gas load produced by this purge is assumed to be a constant Q t . 

Subscript A represents TV Chamber A parameters. Once chamber valves are opened, the overall chamber pressure level 
p A will increase at a rate of 0.0125 Torr/s. The initial chamber pressure is assumed to lie at an insignificant 10' 6 Torr. 

Using the development presented above, the equations used to describe ISIM overpressure in each stage beginning with 
high vacuum conditions is given by 


Vi-%- = Qi-F( P i - Pa ) 

, dt (6) 

= + f {pi ~ Pa) 

at 

For each of the stages identified during the chamber refill process, various physical constraints will produce different 
expressions for gas loads Q and conductance F linking ISIM venting to Chamber A. Note for Chamber A that there is 
no venting term because pumping is disengaged, and the venting term describing ISIM pressure feeds into this volume. 


In the following development, parameters have been chosen for illustrative purposes only, and do not accurately 
represent the ensemble associated with the actual ISIM purge/repressurization application. 



3. PIECEWISE OVERPRESSURE SOLUTIONS 


3.1 Molecular Flow 

The first stage of this process is assumed to occur under molecular flow conditions (free molecule-FM). FM flow is 
governed by Knudsen number Kn, a dimensionless ratio that compares the mean free path length X, a characteristic 
distance between intermolecular collisions, with some characteristic distance d between surfaces. 4 For continuum 
conditions in evacuated vessels, an approximate limit is Kn < 0.01. Theoretically, FM conditions dominate for Kn > 1.0. 
In this application, we will arbitrarily assume FM conditions apply until Kn = 0.01 at the ISIM aperture, where d e g = 

30 cm, ignoring variations that could occur in the transition regime described by 0.01 < Kn < 1. Based on the hard 
sphere intermolecular potential with a molecular diameter of 3.75 A for molecular nitrogen, 4,5 the crossover to 
continuum conditions will occur at X = 3.0 mm, which corresponds to a pressure level of 0.017 Torr. 

For the first hour, repressurization will occur via instrument purges alone, so Q,\ = 0 in (6). Assume the purge rate 
occurs at 1200 L/hr at a supply pressure of two atm, so Qj = 506.6 Torr L/s. If V A = 10000 m 3 while V/ = 40 m 3 : 


at at 


Pa = 


Qi 


(7) 


The chamber repressurization rate under these circumstances is about p A = 5.1 x 10' 5 Torr/s. Substituting into the 
equation for ISIM, we find 


In FM flow, conductance F = A-JrT /2 

t 

Pl(t)=Plfi e 1 


V I ^- = Q I ~F{p I -p A t). 
at 


k is not a function of pressure. The solution for (8) is 


(<1-Pa)i 


l-e 


+ p A t, where (r, q) = 


f V_L Q ± ; 

v F 'Vi 


( 8 ) 


(9) 


v j 

Eq. (9) is meant to describe the passage of the environment around the ISIM aperture from true molecular flow all the 
way to the continuum regime. Even so, the Kn = 0.01 aperture condition is violated a mere 1.5 s after purge initiation! 


3.2 Continuum— Sonic Constraint 

Once the configuration reaches continuum conditions, ISIM venting is constrained by sonic conditions at the AOS 
aperture (Mach number M = 1), and is not affected by chamber conditions as long as the pressure level there remains 
below a certain critical fraction of the ISIM interior pressure created by the purge gases. Assuming isentropic 
conditions, 6 


El 

Pa 


1 + — — ^-M 2 


r 

r - 1 


( 10 ) 


For diatomic gases such as molecular nitrogen at room temperature, the ratio of specific heats y = 7/5. For sonic 
conditions, the maximum chamber pressure that can be supported occurs when M = 1 and p A =0.528 pj . With a 
constant purge flow rate, we may still assume p A = p A t . In Eq. (6), the ISIM venting term is replaced by 6 
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2/, out = mRT = Pi A h RT 
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l2(r-i) 
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Substituting into Eq. (6), 
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dp / 

</f 




The solution becomes 
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where 



( 12 ) 


(13) 


The initial pressure for this stage (subscript 02) is given by the Kn = 0.01 value at the end of the previous FM stage, and 
the At represents time elapsed from then. These conditions hold until 

Pa, 02 + Pa^ = 0-528 p , , where p, = f(At) . (14) 


Calculations performed for these example conditions indicate sonic conditions last for about 6.3 minutes after purge 
initiation, when pj = 0.036 Torr. Beyond this point, rising chamber back pressure influences ISIM internal pressure as 
conditions become subsonic. 


3.3 Subsonic, Instrument Purge Only 

Chamber pressure affects ISIM interior pressure through the vent conductance term in Eq. (6) for subsonic conditions: 6 
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It appears best to solve Eq. (6) numerically. Using superscript n to denote conditions at elapsed time At = nSt : 


Pj(nSt)= p" 1 
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(16) 


where 


p , l(nSt)= p A 03 + p A nSt . (17) 

Initial conditions in this case (subscript 03) are taken from the previous stage at the instant sonic conditions were 
breached. This calculation was performed out to the one hour mark. Using a constant timestep, calculations had to 
iterate in order to produce an acceptably small timestep ( <57 = 0.2 s) so that the last term in Eq. (11) would not begin 
causing the entire solution to become unstable, oscillate, and quickly produce physically unrealistic values. By the end 
of this calculation, chamber pressure has risen to 0.1834 Torr, which is only 3.2e-5 psig below that for ISIM. 



3.4 Subsonic, Instrument Purge & Chamber Valves 

As just mentioned, by the end of that first hour before chamber valves are opened and a new source has been added to 
the system, calculations indicate that overpressure A p = p l — p A has already become a small fraction of either pressure 
alone. This fraction will decrease even further as direct chamber repressurization comes into play, and requires smaller 
timesteps to prevent the instability discussed above. However it turns out that a tremendous simplification can come into 
play when A p/p <0.1. 

In small disturbance theory for s = Ap/p « 1, one may use approximations such as (I + c) a ~ I + a 6 and neglect Of A) 
and higher to greatly simplify Eq. (15): 7 


F sm =F(Ap«p A )*Aj2RT 2A. 
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(18) 


Notice in Eq. (18) that conductance F sm is a function of both p / and p A , but in this incompressible small disturbance 
limit it is independent of y . ; hence it no longer depends on the atomic structure of the gas species under consideration. 

Also, the chamber fill rate has increased considerably from about 5.1 x 10’ 5 Torr/s to p A 4 = 1 Torr/min = 1.25 x 10" 2 
Torr/s, overwhelming the effect of the purge. For small disturbances, one may assume p = p I4 = p 4 4 , which is 
constant. This turns the non-linear ordinary differential equation for p / in Eq. (6) into an algebraic statement that may 
be solved analytically. 


The solution for p I becomes 


V/P = Qi~ A42RT ./paIp/ ~Pa) 
Pa(^)= Pa, 0A+ P A t 


P1 =PA + ^(^P) 2 . 
2 RT Pa { A J 


(19) 


( 20 ) 


Equations 19 & 20 were used to describe the chamber environment up to the point where p A = 760 Torr and 
repressurization was completed. Notice the gas load terms in the numerator of Eq. (20), which represent the strengths of 
sources within ISIM versus that outside of it. This condition must be evaluated to determine if it matches the proper 
physics. It turns out that 506.6 -500Torr L/s is positive (barely), which means the net flow is out of ISIM. Although 
Eq. (20) would still be positive if the external source were higher than the inner one, Eq. (18) would need to be recast in 
terms of p { and A p = p A —pj to reflect the fact that gas would be flowing into ISIM. 

The results described in these sections spanning four flow regimes are depicted together in Fig. 3 below, where net 
overpressure of ISIM above the chamber background is depicted versus time. Early on the overpressure increases with 
time as the ISIM aggregate instrument purge fills that volume with purge gas to higher pressure levels ahead of its 
influence on the chamber, although the maximum overpressure achieved is never quite a sensational level (6.8 x 10" 4 psi 
at 15.5 s after purge initiation). Then the chamber pressure reaches a level where a sonic orifice conditions can no 
longer be maintained, and soon afterwards the pressure differential settles into an inverse time dependence in subsonic 
flow. Once the chamber repressurization valves are opened, the ISIM overpressure amplitude falls below the minimum 
level detectable by ear (about 3 x 10" 9 psig). 10 
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Figure 3. Notional repressurization sequence overpressure spanning four flow regimes. 


4. CREEPING FLOW ANALYSIS 


4.1 Creeping Flow Assumption 

Since we are investigating particulate matter having length scales on the order of microns immersed in low speed flow, 
we expect the Reynolds number based on particle diameter to lie at or below Re = 1, for which creeping flow effects 
should apply. 2 Using room temperature values for molecular nitrogen (dynamic viscosity // of 1.78 x 10" 3 kg/m/s, p = 
1.14 kg/m 3 at 1 atm), 9 


Re =P^< 1 -> (UD) cnt < 


|l.56xl0 5 m 2 /s at 1 atm 

[o.065 m 2 /s at 0.1 834 Torn (see Sect. 3.3). 


( 21 ) 


Equation (23) reveals the condition for creeping flow, which depends on p ^ . If the aperture velocity estimate was 10 
cm/s, the range of validity for creeping flow would include particles below D = 160 microns in diameter at one 
atmosphere, but 65 cm at the initiation of chamber valve opening! Hence, the creeping flow assumption appears to 
cover the ranges of anticipated flow velocities and particle sizes for this application. 


4.2 Net Velocity Calculation 

Once chamber valves are opened, the flow created by this influence is known to stir up particulate matter. Some of these 
particles could potentially settle inside ISIM after falling in through the vertical aperture. Since the reasoning behind 
first applying an aggregate instrument purge is meant to reduce the amount of matter entering that volume later, it was 
decided to convert the overpressure at the aperture to a net velocity by applying an energy conservation statement where 
ISIM acts like a reservoir (velocity U j ~ 0 within the volume for cross-sections much larger than the aperture) and 
potential energy differences in height working against gravity may also be neglected. It is further assumed this statement 
applies at each instant in time quasistatically for this incompressible fluid. Across the aperture (subscript ap), as far as 
these assumptions apply, the statement simplifies to the Bernoulli equation: s 

*P = Pi ~ Pa=^ P u lp ~\p U i ~ 2 pU “P ' (22) 

As this statement is assumed to apply at each instant in time, gas density p may be calculated using the ideal gas law. 
Based on the small disturbance development, the aperture velocity estimate becomes 

U ap = ■JlRTs . (23) 


4.3 Critical Particle Size 

A force balance for lofting a particle with mass m vertically along the y direction with viscous drag in creeping flow 
against gravity (acceleration -g) produces 


my = F S - mg . (24) 

For spherical particles, the viscous drag force is given by the Stokes formula. 2 This formula remains approximately 
valid for a variety of shapes." 


F s =6 7rpR{u ap - y) . 


(25) 


For critical conditions, a particle having an equivalent radius of R cr j t is suspended without vertical motion when 
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P 6 npR c , 


= gT . 


(26) 


In the general solution to Eq. (25), z is identified as an inertial time constant associated with accelerating the particle. 
For an aluminum particle having a diameter of 100 microns, z = 84 milliseconds. Thus, it appears particle motion should 
react rapidly to changes in flow conditions and the quasistatic assumption is justified. 

Equation 26 may be written in terms of R cr j t . For a particle having average mass density p s , the critical radius 
becomes 


Rent — 


f 9 P U ap 
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(27) 


Values of R cr j t and net velocity U ap were plotted versus time from the point where chamber valves were opened in 

Fig. 4 below. A surprising result based on the use of Eq. (25) is that the particle drag force is not a function of gas 
density, unlike the case at higher Reynolds numbers. One consequence of this feature is that larger particles may be 



lofted at low density levels because higher overpressures translate to larger velocities U ap ! Another consequence is that 
the worst case condition occurs at the end of this phase as the chamber just reaches one atmosphere. 
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Figure 4. Maximum suspended particle sizes above ISIM aperture during TV repressurization with chamber valve activated. 


For particle mass density values similar to water and aluminum ( p s = 1.0 & 2.7 g/cm 3 ), 9 predictions indicate the largest 
particle sizes that may be suspended above the ISIM aperture against the chamber valve influence during 
repressurization are 2.0 microns and 1.2 microns, respectively. Since particle fallout distributions are heavily skewed 
towards high concentrations of small elements, these results indicate that it is possible to reject a large fraction of a 
fallout ensemble by number; however, these deflected particles constitute a relatively small fraction of potential area 
concealable by such distributions. 


5. CONCLUDING REMARKS 

An analysis was performed to estimate the ability of an instrument purge to keep particulate matter from entering the 
JWST Integrated Science Instrument Module during chamber repressurization using the facility system. A series of 
models were developed to describe net overpressure across the instrument module’s aperture from molecular flow to 
continuum conditions, with the process also spanning sonic, subsonic compressible, and finally incompressible 
environments. This overpressure was translated to determine aperture velocity for evaluating viscous drag on various 
particle sizes using creeping flow conditions. 


Although the complex set of example conditions used to convey the overall application did not result in a robust ability 
to deny particulate matter entrance to the instrument module by the end of the repressurization process, this situation 
may be remedied by increasing the aggregate purge rate fif possible) or by slowing down the chamber repressurization 
rate. 
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